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Abstract 

The Subgraph Isomorphism problem asks, given a host graph G on n vertices and a 
pattern graph P on k vertices, whether G contains a subgraph isomorphic to P. The re- 
striction of this problem to planar graphs has often been considered. After a sequence of 
improvements, the current best algorithm for planar graphs is a linear time algorithm by 
Dorn (STAGS TO), with complexity 20('=) • 0(n). 

We generalize this result, by giving an algorithm of the same complexity for graphs 
that can be embedded in surfaces of bounded genus. At the same time, we simplify the 
algorithm and analysis. The key to these improvements is the introduction of surface split 
decompositions for bounded genus graphs, which generalize sphere cut decompositions for 
planar graphs. We extend the algorithm for the problem of counting and generating all 
subgraphs isomorphic to P, even for the case where P is disconnected. This answers an 
open question by Eppstein (SODA '95 / JGAA '99). 

1 Introduction 

The Subgraph Isomorphism problem asks, given a host graph G on n vertices and a pattern 
graph P on k vertices, whether G contains a subgraph isomorphic to P. This is a well-studied 
problem that generalizes many other important problems, such as finding cliques, determining 
the girth, finding complete bipartite subgraphs, and finding a Hamilton path or cycle. See 
Eppstein [T3] for a survey on previous results for this problem, and its many applications. 
This problem is NP-complete in general, even for planar graphs (it generalizes Hamilton Path). 
However, in many cases P can be considered to be a small fixed graph. In that case, a trivial 
polynomial time algorithm of complexity n'^^'^^ exists. For general graphs, nothing better is 
known. When restricting G to be planar, this can be improved significantly: Eppstein [131 [H] 
gave a linear time algorithm for Planar Subgraph Isomorphism for any fixed graph P on k 
vertices. This seems best possible. However to judge the practicality of such an algorithm, 
the dependency of the complexity on the value k is also essential. This is where the refined 
viewpoint of parameterized complexity proves to be useful [121 (H]- We choose k to be the 
parameter. An algorithm for this problem is then Fixed Parameter Tractable (FPT) if its 
complexity can be bounded by f{k) ■ 0{n'^), where c is a constant independent of k, and f{k) 
is an arbitrary computable function. 

When viewed as a parameterized problem, the complexity of Eppstein's algorithm |14j is 
20(fciogfc).Q(^^^ hence it is an FPT algorithm. This improved on previous algorithms for Planar 
Subgraph Isomorphism by Plehn and Voigt [20] of complexity 2*^('^'°sfc) . ^O(Vk) ^ ^^^^ Alon et 
al. [2], of complexity 2'^^'^^ ■rpi'^) (which are not FPT algorithms). Finally, Dorn [8] improved 
the previous results and gave an algorithm of complexity 2'-'^^^ ■ 0{n). Eppstein [131 115| also 
considered graphs of bounded genus, which generalize planar graphs. For an integer 5 > 0, 
a graph G has genus at most g if it can be embedded on an orientable surface of genus g, 
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i.e. a sphere with g handles. In |15| . an algorithm for Subgraph Isomorphism in bounded 
genus graphs of complexity 2'^^'^'^°^'^^ • 0(n) is given. In addition, Eppstein [15] considered 
the even more general graph class of apex-minor free graphs, and gave an f(k) ■ 0{n) time 
algorithm, where f{k) is a rapidly growing function of k. We remark that for general graphs, 
FPT algorithms for Subgraph Isomorphism are unlikely |12j . The aforementioned results by 
Eppstein [131 [14} [T5] in fact hold for the more general counting version of the problem, where 
the number of subgraphs of G that are isomorphic to P should be computed. In addition, 
in the case P is connected, he gave an algorithm for listing all of these subgraphs in time 
20(fclogfc) . Q(ji^ _|_ j^-^ . ]^0{i) ^ where m is the number of such subgraphs. 

In this paper, we give an algorithm of complexity 2^^^^ ■ 0{n) for the counting version of 
Subgraph Isomorphism, for the case where G has bounded genus. This generalizes the result for 
planar graphs in [8] and improves the complexity of the bounded genus result in |15j. In addition, 
we give an algorithm that lists all m subgraphs isomorphic to P in in time 

This also holds for the case where P is disconnected, and therefore answers an open question 
from Eppstein |13[ [14] . This is achieved by using a simpler method for counting disconnected 
subgraphs. Our results hold for graphs of bounded non-orient able genus as well. For simplicity, 
we describe the orientable case only, and discuss the non-orient able case in Section [71 In 
Section [7] we also argue that our results apply to finding induced subgraphs. 

There are many examples of problems that can be solved faster on planar graphs and other 
sparse graphs classes such as bounded genus graphs and ff-minor free graphs, in particular 
in the area of FPT algorithms. For instance, for the aforementioned graph classes, many 
parameterized problems can be solved in subexponential time 2^^^^ ■ 0{n), see e.g. [HIZIHII^. 
The theory of bidimensionality [U [5] , easily gives subexponential time FPT algorithms for many 
problems restricted to the aforementioned graph classes. We note that this does however not 
apply to Subgraph Isomorphism, the reason is that bidimensionality applies to problems for 
which the existence of solutions does not change drastically when contracting edges. This is 
certainly not true for Subgraph Isomorphism, apart from a few very special cases for P, such 
as mainly (sets of) paths. 

An essential ingredient for many of these algorithmic results on planar graphs, bounded 
genus graphs and i/-minor free graphs is dynamic programming over tree decompositions and 
branch decompositions. (These are closely related; the minimum width of a tree decomposition 
and a branch decomposition differ by at most a small constant factor.) For problems such as 
finding a minimum independent set or maximum dominating set, straightforward algorithms of 
complexity 2*-^^"'^ • 0(n) exist, where w is the width of the given decomposition, see e.g. |19j . 
For more complex problems such as finding long paths and Subgraph Isomorphism, the best 
known dynamic programming algorithms have a complexity of 2'^("'^°§'") ■0{n). However, when 
restricted to sparse graph classes, this can often be improved to 2^^'^^ ■ 0{n). In the case of 
planar graphs, an essential tool is given by a special kind of branch decompositions, called sphere 
cut decompositions. These were introduced by Seymour and Thomas [22], and their algorithmic 
usefulness was first demonstrated by Dorn et al. |10[lllj. These have subsequently been applied 
many times for constructing fast algorithms. Loosely speaking, a branch decomposition for a 
graph G consists of a labeled tree T, and every edge e G E{T) partitions the edges of G into two 
graphs Gi and G2. In a sphere cut decomposition, for every e € E{T) a simple closed curve in 
the plane exists (a noose), that separates the plane into two regions, one containing Gi and the 
other containing G2. There are many problems that can be solved in time 2<^("''°s"') • 0{n) for 
general graphs on n vertices, when a branch decomposition of width w is given, but that can be 
solved in time 2^^"^^ ■0{n) in the case of sphere cut decompositions [TT]. It is generally believed 
that many algorithms for planar graphs can be extended to graphs of bounded genus. However, 
in the past, making this step has always been a complex technical task. For instance, Dorn et 
al. [9] consider the Hamilton Cycle problem and related problems on graphs of bounded genus, 
and reduce this case to the planar case by cutting the surface a number of times along nooses. 



For the remaining planar case, dynamic programming over sphere cut decompositions is used, 
but this is relatively complex because the previous cuts need to be taken into account. Rue et 
al. [21] proposed a different dynamic programming method for graphs on surfaces: they define 
surface cut decompositions, where two subgraphs Gi and G2 defined by an edge of the branch 
decomposition are separated by a limited number of nooses, which have a limited number of 
common points. 

New techniques and overview of the paper In Section [3] we give a dynamic programming 
algorithm for Subgraph Isomorphism that works for all graphs, when a branch decomposition is 
given. However, in order to give a good bound on its complexity, we restrict to bounded genus 
graphs and introduce a special kind of branch decomposition. One of the main contributions 
of this paper is that we introduce surface split decompositions for graphs of hounded genus, and 
techniques for using these. This is a type of branch decomposition that directly generalizes 
sphere cut decompositions. It allows algorithms that are significantly simpler than previous 
dynamic programming algorithms for bounded genus graphs. In fact, our algorithm and analysis 
is even significantly simpler than that of various previous algorithms for planar graphs (such as 
the previous algorithm for Planar Subgraph Isomorphism [8] ) . Informally, our basic but crucial 
observation is that for surfaces of higher genus, it is irrelevant that the two subgraphs Gi and 
G2 defined by an edge of the branch decomposition can share a complex boundary that consists 
of many nooses; it is only relevant that there are two disjoint (connected) regions Ri and R2 in 
the surface such that Gi is drawn in Ri for i = 1,2. For a precise definition, and an example of 
how this property can be used for giving a good complexity bound, see Section HI In Section [5] 
we give an algorithm for finding surface split decompositions in linear time, by generalizing a 
result by Tamaki [23] and Dorn [8], for finding low width sphere cut decompositions for planar 
graphs of bounded diameter. This is then applied in Section [6] to prove our main algorithmic 
results. We expect that the notion of surface split decompositions, our algorithm for finding 
them, and our technique for bounding the size of corresponding dynamic programming tables 
will be an impetus for the algorithmic research on bounded genus graphs. We believe that 
it will enable the generalization of various existing algorithmic results for planar graphs, and 
that it should allow for the simplification of various known results for bounded genus graphs, 
and more general graphs. In a subsequent paper, we will demonstrate this by applying our 
surface split decomposition techniques to other problems. We remark that algorithms that are 
conceptually simpler are not only convenient for the reader or programmer; often they are also 
faster in practice. Indeed, when restricted to planar graphs, we show in Section [7] that the 
hidden constants in the complexity bound of our algorithm are significantly smaller than those 
in the algorithm by Dorn [8] (although they are still huge). We start by giving definitions in 
Section [2l 

2 Preliminaries 

For basic graph theoretical notations not defined here we refer to [6]. The main graphs that we 
will consider throughout will be simple, but we will construct auxiliary graphs that may have 
parallel edges and loops, i.e. that may be multi-graphs. A walk of length /c in a (multi-)graph 
G is a sequence fo, 61,^1,62, ■ ■ ■ ,Vk where Vi € V{G), and Cj G E{G) is an edge of G that is 
incident with Vi-i and Vi. This is also called a walk from vq to v^. Edges e that are incident 
with u and v will also be denoted by uv if there is no cause for confusion. We will also denote 
walks simply by their vertex sequence vq,vi, . . . ,Vk, which is unambiguous in simple graphs. 
A walk of length is a cycle or k- cycle if vq = Vk, the vertices vq,... ,Vk-i are distinct, and 
all edges are distinct. The distance from u to is the minimum length over all walks from u 
to V. The eccentricity of a vertex u G V^(G) is the maximum distance from u to v, over all 
V G V{G). By d{v) we denote the degree of f G V{G), which is the number of incident edges. 



The operation of subdividing an edge e = uv consists of introducing a new vertex x and edges 
xu and xv, and deleting e. The reverse operation is cahed suppressing the vertex x (which can 
be done when d{x) = 2). An isomorphism between two simple graphs Gi and G2 is a bijective 
function </> : V{Gi) V{G2) such that uv G E{Gi) if and only if (j){u)(l){v) e E{G2). The 
Subgraph Isomorphism (counting) problem takes as input a simple (host) graph G and a simple 
(pattern) graph P. The objective is to compute the number of subgraphs G' of G that are 
isomorphic to P. Such a subgraph G' is called a P-isomorph, or isomorph if the graph P is 
clear. Observe that the case where both G and P are non-simple is easily reduced to the simple 
case by subdividing every edge. 

Branch Decompositions The set of leaves of a tree T is denoted by L{T). A branch 
decomposition of a graph G is a tuple (T, //) consisting of a ternary tree T (i.e. all non-leaves 
have degree 3), and a bijection ^ : -i^(T) E{G). For a subset 5 C L{T), we will use /u(«S') to 
denote the set of images. Every edge ct G E{T) defines a middle set mid{eT) ^ ^(G) in the 
following way: let Ti C V{T) and r2 C y(T) denote the vertex sets of the two tree components 
oiT — eT- The edge sets ii{L{T)r\Ti) and fi{L{T)r\T2) partition the edges of G. By Gi and G2 
we denote the respective subgraphs of G induced by these edge sets. Now mid(eT) is defined 
as V{Gi) n V{G2)- Hence this is the set of vertices that are both incident with edges belonging 
to Ti and edges belonging to T2. The width of (T, fi) is defined as 

width(r, /i)= max |mid(eT)|- 

eT&E{T) 

A rooted branch decomposition is a tuple (T, fj,) where T is a ternary tree, and a root r € L{T) is 
identified. In this case, /i is a bijection from L{T)\{r} to E[G). A rooted branch decomposition 
can easily be obtained from a branch decomposition by subdividing an arbitrary edge with a new 
vertex x, and adding the vertex r and edge rx. In the case of a rooted branch decomposition, 
for e € E{T), by Tg C y(T) we denote the set of vertices of the component of T — e that does 
not contain r. Similar to above, Tg defines a subgraph of G, which is denoted by Gg. Since T 
is ternary, every edge e G ^(P) for which Tg is non-trivial (i.e. not a single vertex) has two 
children; these are the other two edges of T that are incident with the end vertex of e in Tg. 
Observe that if e,. G ^(P) is the edge incident with r, then Gg^ = G. 

Graphs embedded on surfaces For an introduction to graphs embedded on surfaces we 
refer to |18j . Formally, a surface is a connected compact 2-manifold without boundary. For 
every integer g > 0, Let Sg denote a surface that is obtained by adding g handles to a sphere. 
Hence §g is an orientable surface of genus g, and in fact, for every orientable surface there 
is exactly one value of g such that it is homeomorphic to Sg for [18^ p. 82]. For ease of 
presentation, all surfaces that we consider will be orientable. In Section [7] we discuss how to 
extend our proofs to the non-orientable case. A region of a surface is a connected open set. The 
boundary of a region R consists of all the points that lie in the closure of R but not in R itself. 
For a simple curve C Q Sg, all the points in G that are not the end points are called interior 
points. 

An embedding of a graph G into the surface Eg consists of an injective mapping of the 
vertices of G to points ipiv) in §g, and a mapping of edges e = uv of G to a simple curve 
ip{e) in Sg with end points ip{u) and ip{v), such that i{j{e) shares no interior points with ^(e') 
or 'ip{v) for any e' G E{G)\{e} or -y G V{G). (That is, edges do not cross and only their end 
points overlap with vertices.) To simplify terminology and notation, if tp is an embedding of a 
graph G, then the images ipiv) and ip{e) (which are subsets of Sg) for v G V{G) and e G E{G) 
will also be called vertices and edges of G, respectively. Let X C Sg he the union of all vertices 
and edges of an embedded graph G. The faces of G are the maximal regions of Sg\X. An 
embedding is a 2-cell embedding if every face is homeomorphic to an open disc. For a graph G 
with n vertices and m edges, which is 2-cell embedded in E>g with / faces, by Euler's formula it 



holds that n — m + f = 2 — 2g |18| p. 85]. The boundary of every face F oi G defines a closed 
walk 1^ in G in a straightforward way, which will be called the facial walk for F. If the length 
of W is fc, then F will be called a k-face. 

Given a 2-cell embedding ^ of a graph G in §g, we can define the following two related 
(multi-) graphs, and their 2-cell embeddings in §g. The dual graph G* of G has one vertex vp 
for every face F of G, and one edge e* for every edge e € E{G), with end vertices vp and vp', 
where F and F' are the faces incident with e. Note that G* may have parallel edges and loops. 
e* is called the dual edge of e and e the primal edge of e*. A straightforward 2-cell embedding 
ip* of G* exists such that vp lies in the face F, and such that for e € E{G) and e' G E{G*), 
ip{e) overlaps with ^*(e') if and only if e' is the dual edge of e. 

The radial graph TZq of G (which is also called the vertex-face incidence graph) is obtained by 
starting with the vertex set V{G) (the original vertices), and adding a vertex vp for every face F 
of G (the face vertices). For every face F in G with corresponding facial walk uq, tti, . . . , uq, 
add k edges uqVp, uivp, . . . , u/^^ivp. Note that T^g ™ay have parallel edges again. The original 
vertices and face vertices of TZq define a bipartition of TZq- A straightforward 2-cell embedding 
tp* of TZq exists such that for every face F of G it holds that ip*{vp) € F, for every v G V{G) 
it holds that ip*{v) = ip{v), and for every edge e = uvp of IZq, interior points of ip*{e) lie in 
the face F. 

A combinatorial embedding vr of a graph G consists of a cyclic order tt^ on the incident 
edges for every vertex v G V{G). Prom an embedding of a graph we obtain the corresponding 
combinatorial embedding by considering the clockwise order of edges around every vertex. A 
map is a connected graph G together with a combinatorial embedding vr. The set of all facial 
walks and therefore the number of faces of a map can easily be deduced (without constructing 
the embedding), so the number of faces / of a map is well-defined. The genus g of a map G, vr is 
the solution ton — m + / = 2 — 2g, where n = \V{G)\, m = |£'(G)| and / is the number of faces 
of G, IT. The genus of a connected graph G is the minimum g such that G admits a combinatorial 
embedding of genus g. Given a map G, tt of genus g, a corresponding 2-cell embedding of G 
in Sg exists. Our algorithms only work with combinatorial embeddings; it is not necessary to 
generate an embedding in Eg. However for the correctness proofs it is convenient to argue with 
corresponding 2-cell embeddings of G in Sg. 

3 A dynamic program for counting colorful isomorphs 

In this section we give a dynamic programming algorithm for the following generalization of 
Subgraph Isomorphism. An instance of the Colorful Subgraph Isomorphism problem consists of 
a colored host graph G, and a pattern graph P. The coloring of G is a function a : V{G) — ?> G, 
with G = {1, . . . ,q}. This encodes a partition of V{G) into q sets. We remark that this is 
not required to be a proper vertex coloring, so adjacent vertices may receive the same color. A 
subgraph G' of G is called colorful if for every color x ^ C, G' contains a vertex v of color x. 
(Note that G' may have more than q vertices.) The objective is to count the number of colorful 
P-isomorphs of G. We now present an algorithm for this problem. (When q = 1, this is the 
original counting problem.) 

Dynamic programming table Let (T, ^u) be a rooted branch decomposition of G. For 
every edge e G E{T), we will form a dynamic programming table %■ Informally, this table 
will store information about all possible subgraphs of the graph Gg, on at most k vertices. 
Firstly, we distinguish between non-isomorphic subgraphs. Furthermore, subgraphs of Gg that 
are isomorphic but intersect differently with the 'boundary' mid(e) of Gg are also considered 
distinct. Finally, we keep track of the set of colors that appear in these subgraphs. Two 
subgraphs of Gg are only considered equivalent if they match in all three regards. In that case, 
there will be a single entry in the table that represents both subgraphs. We now define this 




Figure 1: An illustration of the dynamic programming update step: An edge e € E{T) with 
children f,g(z E[T) defines the subgraph Ge of G. A subgraph Gg of Ge is shown in bold. A 
table entry {He, 'ye,Ae,r]) € Te that corresponds to Gg is obtained by combining the compatible 
entries {Hf,'yf,Af,T]f) G 7/ and {Hg,jg,Ag,T]g) G Tg. Vertex colors are ignored in the figure. 

formally. The definitions and dynamic programming update step are illustrated in Figured] 

Let -ff be a graph, and let 7 be a mapping from mid(e) to V{H) U {nil}, which is injective 
on V{H). To be precise, every vertex of V{H) occurs at most once as a 7-image, but multiple 
vertices may be mapped to nil. Furthermore, let A C {1, . . . For such a tuple {H,^,A), a 
subgraph G' of Gg is called an {H,^, A) -subgraph if the following two properties hold. 

• There is an isomorphism <j) : V{G') — )• V{H) with 7(f) = (j){v) for all v G mid(e) fl V{G'), 
and 7(7;) = nil for all v G mid(e)\y(G'). 

• For all colors x G G, x G ^ if and only if G' contains a vertex of color x. 

For e G E(T), the dynamic programming table Te will now contain entries {H,^, A,r]), where 
H, 7 and A are as defined above, and 7/ is a non-negative integer. The idea is that such a 
table entry indicates that Gg contains exactly r] non-equivalent {H, 7, A)-subgraphs. Two table 
entries (-ffi, 71, ^1, ryi) and (-^25 72, f?2) are equivalent if the following properties hold. 

• There is an isomorphism (j) : V{Hi) — ?> V{H2) such that for all v G mid(e), either 'yi{v) = 
72 (^) = nil, or (j){'yi{v)) = 72 (t') holds. 

• Ai = A2. 

Observe that the above definition satisfies the following property: a subgraph G' of Gg is 
both a (ffi, 7i, Ai)-subgraph and a (//2 5 72, ^2)-subgraph if and only if (ffi , 71 , , 7/1 ) and 
(i?2, 72, ^2, are equivalent. Therefore, when two entries (i^i , 71 , Ai , 7/1 ) and (-^21 72, ^2, ''?2) 
are equivalent, we can merge them by replacing them with the single entry 71, ^1, 771 +7/2). 
We say that the table Te is k-correct if 

• for every tuple (H, 7, ^), T contains an entry {H, 7, A, rf) if and only if Gg contains exactly 
7/ > 1 graphs G' with |y(G')| < k that are (ff, 7, A)-subgraphs, and 

• Te contains no pairs of equivalent entries. 



Dynamic programming update step Let e G E{T) be an edge with children / and g. 
We will now show how a fe-correct table Te for e can be obtained from fe-correct tables 7/ and 
Tg for / and g, respectively. Let and G'g be subgraphs of G/ and Gg, respectively. They 
are called compatible if for every v G mid(/) n m[d{g) it holds that v G V{G'j) if and only 

if v G V{G'g). Note that the function V(G'^,G^) := (v{G'f) U ^(G;,), U -E(G;^)) is a 

bijection between pairs of compatible subgraphs of Gj and Gg and subgraphs of Gg. Similarly, 
we define two entries {Hf,jf,Af,r]f) G 7/ and {Hg,^g, Ag,r]g) G 7^ to be compatible if for all 
f G mid(/) n mid(5), it holds that 7/(v) = nil if and only if ^g{v) = nil. We now show how to 
combine two such compatible entries into an entry (iJg, 7e) ^ej ^e) that satisfies the following 
property: if G'j is an {Hf, 7j, Ay)-subgraph and G'g is an {Hg,^g, Ag)-subgraph, then ip{G'jr, G'g) 
is an (ffg, 7e, j4e)-subgraph. 

• For all V G mid(/) PI ui[d{g) with 7/(v) ^ nil (and thus 7g(u) 7^ nil): identify the vertex 
Jf{v) of i/j with the vertex jg{v) of ffg, and call the new vertex i^{v). This gives the 
graph Hg. 

• For all V G mid(e): 

1. If 1) G mid(/)\mid((7): set Jeiv) = ^f{v)- 

2. If t; G mid((7)\mid(/): set 76(1^) = lg{v)- 

3. If G mid(5') nmid(/): set 7e(t') = i^('t^)- 

By definition of mid(e), this covers all cases and thus defines the function 7e. 

• Set Ae = AfVJ Ag. 

• Set r/e := ?// • 

The above definitions and observations show that if there are r/j (//j, 7^, j4j)-subgraphs in G/, 
and rjg (//g, 7^, j4g)-subgraphs in Gg, then there are r/j • rjg (ffg, 7e, 74e)-subgraphs in Ge that are 
the result of combining graphs of the former two types in every possible combination. However, 
there may also be (//g) 7e; ^e)-subgraphs of Ge that are the result of combining different types 
of subgraphs of G/ and Gg. Therefore, the entire dynamic programming update step requires 
the merging of entries as well. We will now present this in detail. 

Assuming that we have /c-correct tables T/ and Tg for / and g respectively, we construct 
a /c-correct table Te for e as follows: We start with an empty table T- Then we consider 
every pair of compatible entries from 7/ and Tg, and combine them as described above. For 
every combination, this yields a possible entry (H,^, A,r]) for T- In case that H has more 
than k vertices, we ignore this possible entry. Otherwise, we check whether T already contains 
an equivalent entry {H' ,'y' , A' ,r]'). If so, we merge the two entries. If not, we add the entry 
(T, -y,A,r]) to the table T- Then we continue with the next pair of compatible entries from 7/ 
and Tg- 

Lemma 1 Let (T, fi) be a rooted branch decomposition for a colored graph G, and k be an 
integer. Let e G E{T) be an edge with children f and g, for which k-correct tables Tf and Tg are 
given. Then the table T that is constructed with the above dynamic programming update step 
is k-correct for e. The construction takes time ■ f{k) ■ k'^^^\ where f{k) is the complexity of 
deciding whether two entries are equivalent, and X is an upper bound on the size of a k-correct 
table. 

Proof: We argue that the construction of T is correct. Clearly, T contains no pair of equiv- 
alent entries. First we show that for every resulting entry (Hej'je, Ae,r]e) G T, there are at 
least r]e (-ffe, 76, ^e)-subgraphs in Gg. Since 7/ is A;-correct, every entry {I{f,^f,Af,rjf) G 7/ 



corresponds to a set of r]f (//j, 7j, ^j)-subgraphs of Gf. A similar statement holds for Tg. 
Consider compatible entries {Hf,^f,Af,rif) € 7/ and {Hg,^g, Ag^rjg) € Tg. The bijection 
and the construction of (i^e, 7e) ^e, ^e) show that there are r]f • rjg subgraphs of Ge that 
are (ffg, Te, ^e)-subgraphs, that are obtained by combining every (ffj, 7j, Aj)-subgraph with 
every {Hg,jg, Ag)-subgraph. When merging a new entry {He,'ye, Ae,r]e) with an existing en- 
try (i/g, 7g, r^g), the corresponding two sets of subgraphs of Ge are disjoint, because they 
correspond to different types of subgraphs oi Gf and Gg. This shows that for every entry 
(Hej'je, Ae,r]e) € Te, 7?e IS a lower bound for the number of (i?e, 7e) ^e)-subgraphs in Ge- 

We conclude the correctness proof by showing that every graph is counted this way. Consider 
a subgraph Gg of Ge on at most k vertices. Its restrictions to G/ and Gg respectively give a 
compatible pair G'j and G'g, both with at most k vertices. Since 7/ and Tg /c-correct, there 
are unique entries {Hf,^f,Af,rif) E 7/ and {Hg,jg,Ag,r]g) € Tg such that G'^ is one of the 
r]f (//j, 7j, j4j)-subgraphs of Gf, and a similar statement holds for G'g. This shows that the 
subgraph Gg of Gg is counted by the dynamic programming update step; it is one of the rjf ■ rjg 
resulting subgraphs. 

Finally, we consider the complexity. Both 7/ and Tg contain at most X entries, and the 
same holds for T, at every time during its construction. Therefore, there are at most X'^ pairs 
of entries from 7/ and Tg that should be considered. Testing whether they are compatible, and 
if so, combining them can be done in polynomial time k'-^^^\ Subsequently, there are at most 
X entries of T that should be tested for equivalence, which takes times f{k) ■ X. This proves 
the stated complexitjQ. □ 

This dynamic programming update step is the core of the following dynamic programming 
algorithm: 

• First, for every edge e S E{T) that has no children, Gg consists of a single edge, so it 
is trivial to compute a fc-correct table T- To be precise, assuming that mid(e) ^ 0, the 
five entries of T correspond to the empty graph, a graph consisting of a single edge, two 
graphs consisting of a single vertex (with different 7), and one graph on two vertices, 
without edges. 

• For every edge e G E{T) with two children / and g, we compute a /c-correct table T 
using the dynamic programming update step, after /c-correct tables for 7/ and Tg have 
been computed. 

• After computing /c-correct tables for all edges of T, we inspect the table T^ where is 
the root edge of {T,fi). Since Gg^ = G, mid(er.) = 0. Therefore, T^ contains at most one 
entry (H, 7, A, rj) such that H is isomorphic to P and A = {1, . . . , q}. If there is such an 
entry {H,'y, A,t]), we return 1], and otherwise we return 0. 

The correctness of the above procedure follows from the definitions. Combined with Lemma [1] 
this gives the following theorem. Note that a branch decomposition (T, /i) of a graph on m 
edges has \E{T)\ G 0(m). 

Theorem 2 Let (T, n) be a rooted branch decomposition of a colored graph G with m edges. In 
time X^-f{k)-kO^^^-0{m), it can be computed how many colorful subgraphs of G are isomorphic 
to a given graph P on k vertices, where f{k) is the complexity of deciding whether two entries 
are equivalent, and X is an upper bound on the size of a k-correct table. 

The above theorem applies to general graphs, for which an appropriate bound for X can be 
given. However, to obtain the desired complexity, we introduce surface split decompositions in 



^We remark that there is room for improvements in the the complexity (constant factors) here: when the 
tables are ordered, not all entry combinations have to be tested for compatibility, and not all entries have to be 
considered for testing equivalence. See e.g. [19] for examples of such improvements. For simplicity we ignore this 
in our analysis. 



the next section. In the case that an embedding vr of G of bounded genus is given, and (T, fx) 
is a surface split decomposition, we give a good upper bound for X. 

4 Surface split decompositions and a bound for the table size 

In this section, we define surface split decompositions for graphs embedded in E>g. These gen- 
eralize the notion of sphere cut decompositions, which are defined for graphs embedded on the 
sphere So (i-e- plane graphs). Given an embedding of G in the sphere §0) a noose is a simple 
closed curve in Sq that intersects G only in its vertices. Hence by the Jordan Curve Theorem, 
'§>o\N consists of exactly two connected components, the regions i?i and i?2. A branch decom- 
position (T, /i) of G is a sphere cut decomposition [llj if for every ex € E{T) with corresponding 
subgraphs Gi and G2 of G, there exists a noose A'' that divides So into Ri and i?2, such that 
for i = 1,2, Gi lies in U Ri, and intersects G exactly in mid(e). Our simple but essential 
observation is that the fact that the two regions Ri and R2 are separated by a simple closed 
curve A^ is actually irrelevant; this is merely an artefact of restricting to the sphere. We will 
show that the appropriate way of generalizing sphere cut decompositions to surfaces of higher 
genus is surprisingly simple, as follows. 

Definition 3 Let G be a graph embedded onSg. A branch decomposition (T, /i) of G is called 
a surface split decomposition if for every e € E{T) and corresponding subgraphs Gi and G2 of 
G, there are disjoint regions Ri C Sg and R2 C Sg such that for i = 1,2, all vertices and edges 
of Gi are drawn entirely in the closure of Ri . 

Observe that this definition implies that all vertices in mid(e) lie on the boundary of the 
closure of Ri and on the boundary of the closure of i?2- Note also that w.l.o.g. we may assume 
that all other vertices of Gi and internal points of edges of Gi lie entirely in Ri itself (not just 
in the closure). We stress that it is crucial that Ri and R2 are connected open sets. If not, 
then firstly the above definition is not a generalization of sphere cut decompositions, but more 
importantly, the proof of Lemma[7]below fails. It is essential that there is no single point x G i?j 
such that Ri\{x} consists of multiple connected components (for i = 1,2). Finally, we remark 
that even if the regions i?i and R2 are chosen such that their boundaries are the same, this 
boundary is not necessarily a simple curve if 5 > 1. This is illustrated in the figures below. 
The definition easily extends to maps G, vr of genus g: (T, fi) is a surface split decomposition of 
G, vr if it is a surface split decomposition for any embedding of G in Sg that corresponds to the 
combinatorial embedding vr. 

Let (r, fi) be a surface split decomposition for a map G, vr. In this section, we firstly give a 
bound on the size of a /c-correct table, using the fact that (T, /i) is a surface split decomposition. 
Secondly, we show that equivalent pairs can be recognized in polynomial time. 

For the bound on the table size, we will use a bound on the number of different graphs on 
n vertices, embedded in a surface of genus g. To be precise, we will consider simple, connected 
graphs G, that come with a combinatorial embedding vr. Such a pair G, vr is called a simple 
map. In addition, a tuple {u,v) of vertices is given such that uv € E{G). This is the root, in 
Figure [2] indicated with an arc. Such a combination G,tt, {u,v) is called a simple rooted map. 
Two simple rooted maps G, tt, (u, v) and G' , vr', (u', v') are equivalent if there is an isomorphism 
(j) : V{G) — )■ V{G') with (f>{u) = u' , (j){v) = v' , and that preserves the clockwise order of edges 
around every vertex, or reverses all of these orders. In other words, (f> should map faces (facial 
walks) of G to faces of G' . In case an edge labeling is given for both simple rooted maps, we 
also require that (p maps edges to edges of the same label. To bound the number of simple 
rooted maps, we apply a result by Bender and Canfield [3^1 . 



^Bender and Canfield included multi-graphs (i.e. rooted maps) in their bound, but clearly this yields an upper 
bound for the number of simple rooted maps. 



(a) Ge and region R\. (b) Vertex u and labeling of mid(e) (c) H and 7 from an entry 

a subgraph is shown m bold. ^jy ^ ^ 




(d) Embedding H and adding edges, (e) A rooted labeled map J?', A and function p 

that encode H and 7. 

Figure 2: An illustration of the proof of Lemma [71 The rectangle represents a torus; edges that 
leave the rectangle on the left continue at the same height on the right, and a similar statement 
holds for edges leaving the top and bottom of the rectangle. 



Theorem 4 ([3]) There are 0{m^^^ i)/2]^2'") equivalence classes of simple rooted maps of 
genus g on m edges. 

Using an elementary argument based on Euler's formula, one can prove that the following 
bound holds in the case of simple graphs. 

Proposition 5 Let G he a simple graph of genus g, with n vertices and m edges. Then m < 
3n — 6 + Qg. 

Combining these bounds yields the following rough upper bound in terms of the number of 
vertices. 

Corollary 6 Let g be a constant. There are 2^^"^^ simple rooted maps of genus at most g on n 
vertices. 

Proof: When considering g to be constant. Theorem H] and Proposition [5] show that the number 
of simple rooted maps on n vertices of genus g can be bounded by 

O ((3n + 65)^(f^^)/2 . i2^n . 12%^ c n^(5) . 2OW c 2^("). 

Clearly, this is then also an upper bound for the number of simple rooted maps of genus at 
most g. □ 

Lemma 7 Let g be a constant. Let G, vr be a simple map of genus at most g, for which a 
coloring with q colors and a surface split decomposition (T, /i) of width w are given. Let k be 
an integer. For e G E{T), let Te be a k-correct table. Then \Te\ e 2^ ■ 2"" ■ 2^^^\ 

Proof: For all entries (-ff, 7, A^rf) S 7^ where H has at most k vertices, we will encode Li and 
7 by a simple rooted map H',Tr', {u,v) on at most k + 1 vertices, together with a 0/1-labeling 
A of a subset of the edges of H' , and a 0/1-labeling p of the vertices in mid(e). This means 
that any two non-equivalent entries (i^i , 71 , >li , 771 ) and (if2, 72, ^2, either have Ai ^ A2, or 



yield a different labeling p, or yield non-equivalent rooted maps. The encoding is illustrated in 
Figure El 

We use the following auxiliary graph. Since (T, ii) is a surface split decomposition, there 
are disjoint regions i?i and R2 in Sg, such that Ge lies in the closure of and mid(e) lies 
on the boundary of both Ri and R2. Thus we can extend Ge by drawing a vertex u in R2, 
and drawing edges in the closure of R2 from u to every vertex in mid(e), while maintaining 
an embedding in S^. Number the vertices of mid(e) fi, . . . ,vt, corresponding to the clockwise 
order of edges around u E|. 

Now we show how to construct the encoding H',tt', p, A for an entry {H,j, A,rj) G Te- 
Firstly, for all vertices v S mid(e), set p{v) = if and only if 7(f) = nil. The graph H' is 
constructed as follows. Since Te is A:-correct, H corresponds to a subgraph of Ge, so H can 
also be drawn in the closure of Ri, such that all vertices that are 7-images are drawn on the 
boundary of Start with such an embedding. Next, add the vertex n, and edges uvi for every 
Vi G mid(e) with 'y{vi) ^ nil. Draw these as described in the previous paragraph. This yields 
a simple graph embedded in S^, on at most k + 1 vertices. However, it may not be connected 
(if H has components that contain no 7- image). Add edges between different components until 
the graph is connected. This yields H' . Clearly, drawing these new edges can be done while 
maintaining an embedding. Hence H' is embedded in E>g, and the corresponding combinatorial 
embedding vr' has genus at most g. To obtain a rooted map, choose i to be the lowest index 
such that p{vi) = 1. Then the tuple {u,Vi) is chosen as the root of H' ,tt' . That is, v = Vi. A 
bridge of a connected graph G is an edge e G E{G) such that G — e is disconnected. For all 
bridges e G E{H'), we set A(e) = 1 if e G E{H), and A(e) = otherwise. 

To argue that this is an encoding, we show that the pair H, 7 can be reconstructed from 
H' , vr', [u, v),p. A, without using information about vertex labels of H' (since we apply a bound 
on equivalence classes of maps, or in other words, unlabeled maps). Let i be the lowest index 
such that (t){vi) = 1, and let {u,v) be the root of H' . Then ^{vi) = v. Let i' be the next index 
such that (j){vi') = 1, and let {u,v') the the next edge incident with u, in clockwise order. Then 
l{vi') = v' . Continuing this way, and setting j{vj) = nil for all vj with p{vj) = 0, we obtain 
7. Next, delete u and every bridge e with A(e) = 0. This way we obtain the graph H. This 
yields the following conclusion: Suppose two entries (ifi, 71, y4i, 771) and {H2, '^2, ^2,112) yield 
encodings H[,'n'^,{ui,vi), pi,\i and H'2,tt'2,{u2,V2), P2,^2 respectively, such that pi = p2 and 
such that the edge labeled simple rooted maps H[,Tr[, {ui,vi) and H2, 1^2, {u2,V2) are equivalent. 
(Recall that this requires that the isomorphism maps edges of to edges with the same label.) 
Suppose that in addition Ai = A2. Then (i?i, 71, ^1, 771) and (-^2; 72, ^2, '72) are equivalent. 

Since a fc-correct table contains no pairs of equivalent entries, we can now give an upper 
bound for its size. There are at most 2^ possibilities for p, at most 2'^ possibilities for A, at 
most 2'-'^^^ simple rooted maps on at most A; + 1 vertices of genus at most g (Corollary [6]), and 
at most 2*^ possible labelings A (since a graph on /c + 1 vertices contains at most k bridges.) 
Therefore, the number of entries in a A;-correct table is bounded by 2« • 2^" • 2<^('=). □ 

Now we show how to test whether two table entries are equivalent. This is based on the 
algorithm by Miller |17j . for testing whether two graphs of bounded genus are isomorphic. 
See [16] for an alternative algorithm with the same complexity. 

Theorem 8 ([Tf]) In time rP^9\ it can he tested whether two graphs on n vertices of genus g 
are isomorphic. 



''We remark that if R2 is not iiomeomorphic to a disc, then there are multiple ways of drawing these edges, 
resulting in different clockwise orders. Some will yield a better bound than others, but for proving the stated 
upper bound that is not important; it is only important that we choose one way of drawing these edges and use 
this for every entry of Te.- We remark also that at this point we require that R2 is not only connected, but also 
open, in order to draw the edges. 



Proposition 9 Let (T, fx) be a branch decomposition of a map G, vr of genus g, and let Te be 
a k- correct dynamic programming table for e e E{T). In time k^^^) ^ 

it can be decided whether 

two entries (i^i , 71 , Ai , 771 ) and (-f^2, 72, ^2, ^2) of Te are equivalent. 
Proof: We assume that 

• Hi and H2 have the same number of vertices, 

• for every v G mid(e), 71(f) € ^(-^1) if and only if 72('y) € V{H2), 

• for every v G mid(e) with 71 (v) 7^ nil, 71(f) and 72(f) have the same degree, and 

• Ai = A2. 

Otherwise, we may immediately conclude that the entries are not equivalent. We shall make 
small modifications to both H\ and H2, to ensure that any isomorphism (f) of the resulting 
graphs H[ and H2 maps 71(f) to 72(f) for every v £ mid(e). Then it suffices to test whether 
H[ and H2 are isomorphic. 

Number the vertices f G mid(e) with 71(f) ^ nil as fi, . . . ,fj, in an arbitrary order. We 
modify Hi by adding k vertices, all joined by an edge to 71 (fi), and adding 2k vertices joined 
to 7i(f2), etc. This ensures that in the resulting graph H[, for all i, the vertex 7i(fi) has a 
unique degree. Apply the same modification to H2 to obtain H2. Since the table Te is /c-correct, 
both Hi and H2 are isomorphic to subgraphs of G, and therefore have genus at most g. Adding 
degree 1 vertices does not increase the genus, so H[ and H2 have genus at most g as well. 

Because of the unique degrees, if there exists an isomorphism (j) : V{H[) V{H2), then 
for all i, </'(7i(fi)) = 72(fi)- It follows that H[ and H2 are isomorphic if and only if {Hi,-yi) 
and (-^2,72) are equivalent. The former can be decided in time x^^^^ (Theorem [8]) , where 
X = \V{H[)\. Since x G 0{k^), this shows that we can test equivalence in time k'^'^^d) = A:'^(5').n 

Combining Theorem[2]with Lemma[7]and Proposition[9]gives the following theorem. Here we 
use that a graph on n vertices of genus g has at most 0{n+g) edges, and that k^^^^ ■2'-'^^^ C 2'-'^^\ 

Theorem 10 Let g be a constant. Let G,Tr be a simple map on n vertices of genus at most g, 
for which a coloring with q colors and a surface split decomposition (T, fi) of width w are given. 
Let P be a graph on k vertices. In time 8'' ■ S'" ■ 2'^^'^) • 0{n), it can be computed how many 
colorful subgraphs of G are isomorphic to P. 



5 Constructing surface split decompositions 

Tamaki [23] gave a linear time algorithm for constructing a branch decomposition of width 
2d+\ of a graph G embedded in the sphere, when a vertex r G V{G) of eccentricity d is giver0. 
Dorn gave a different presentation of the construction and showed that it yields in fact a 
sphere cut decomposition. We show how the construction of [231 Ej can be extended to surfaces 
of arbitrary genus, and prove the following theorem. 

Theorem 11 Let G, vr be a map of genus g, for which a vertex r G V{G) of eccentricity d is 
given. In linear time, we can construct a surface split decomposition (T, fj.) of G of width at 
most (25 + l)(4d + 3)/2. 

Note that for g = 0, Theorem 1111 implies the aforementioned result from p3] and |8] since 
[(4(i+ 3) /2j = 2d+l. The remainder of this section will be devoted to the proof of Theorem 111! 
Before we construct the surface split decomposition {T,fi), we will construct a number of 
auxiliary graphs. These are illustrated in Figure [3l 



^Tamaki formulates the result differently, such that r corresponds to the outer face of a graph embedded in 
the plane. 



The graph G' He and G' 




An oriented BFS tree Tg of TZc A spanning tree T* of the dual of TZc 




A partition of V{T*) into Ti and Two subgraphs of G' and mid(eT) 
and the corresponding regions. 



Figure 3: An illustration of the graphs used for the proof of Theorem II 11 for the case of 5 = 1. 
The rectangle represents a torus; edges that leave the rectangle on the left continue at the same 
height on the right, and a similar statement holds for edges leaving the top and bottom of the 
rectangle. 



The objective of the first stage of the construction is to construct a tree T*, such that for 
every e € E{G), there is a unique vertex of T* associated with e. This can be thought of as the 
function fi, from a subset of V{T*) to E{G). We will show in Lemmas 1121 and 1141 below that 
this T* is already 'almost' the desired surface split decomposition. However, T* will not yet 
be ternary (though it will have maximum degree 3), and the vertices that are associated with 
edges of G are not necessarily leaves. This is subsequently addressed in the second stage of the 
construction. 

First stage of the construction First we modify the graph G as follows: for every edge 
e = uu, add two extra parallel edges between u and one on either side of e (while maintaining 
a combinatorial embedding). This ensures that all original edges of G are incident with two 
2-faces, and that every vertex has degree at least 3. Denote the resulting embedded graph by 
G' . The edges in E{G') PI E{G) are called original edges, and edges in E{G')\E{G) are called 
new edges. Observe that this operation does not change distances in the graph. (In Figure O 
for the sake of readability we have omitted this first step, and assume that G = G' . This is no 
problem since this doubling step is only used to guarantee degree upper bounds for the resulting 
tree T*.) 

Now construct TZq', the radial graph of G' . Let T5 be a BFS spanning tree of T^g'j rooted 
at a vertex r e V{TIg') n V{G) of eccentricity d. Choose T* to be a spanning tree of the dual 
graph of TZg', such that for all edges e* G E{T*), the corresponding primal edge e is not in Ts- 
Note that it is possible to construct such a spanning tree, since Ts contains no cycles. 

If the genus g of the given map is zero, then at this point every edge of TZq' appears either 
in Ts, or its dual edge appears in T* . However, if g > then there are some edges for which 
neither property holds. Add all of these edges to Ts, to obtain the graph Tg . This completes 
the first stage of the construction of a surface split decomposition. For the proofs below we 
switch to a topological viewpoint, and assume that all of the given graphs are embedded in Sg 
with a 2-cell embedding, as discussed in Section [2j 

We remark that there is a trivial bijection between the faces of TZq' and the edges of G' , 
since every face of TZc contains exactly one edge of G' . Secondly, there is a trivial bijection 
between the vertices of T* and the faces of TZq', since T* is a subgraph of the dual graph of 
TZg'- Combining these bijections yields a bijection from V{T*) to E{G'). In the remainder of 
this section, we refer to these bijections when speaking of 'the set of edges of G' or faces of TZc 
that corresponds to a subset of V{T*y, or 'the set of edges of G' that corresponds to a subset 
of the faces of T^g''- The subgraph of G' that corresponds to a set Ti C V{T*) is the subgraph 
of G' induced by the edges that correspond to Ti . 

For a set T of faces of TZc, we define the open set R{T) C Sg as follows: Informally, R{J^) 
is obtained by starting with T and 'connecting' all pairs of faces that share an edge. To be 
precise, to obtain R{J-) we first delete all edges and vertices of TZc for which all incident faces 
are in T. Call the resulting graph TUq,. Now R{J~) is the union of faces of TZ'q, that contain 
points of faces in T. Observe that if the set T corresponds to a connected subgraph of the dual 
graph of T^G') then R{J^) is connected, and therefore a region. For a set Ti C V{T*), we denote 
by R(Ti) the open set R{J^), where is the set of faces of TZq' that corresponds to Ti. The 
following property now follows from the definitions and observations above. 

Lemma 12 For an edge ex € E[T*), let Ti and T2 be the vertex sets of the two components 
ofT* — ct, and let G'l and G'2 he the subgraphs of G' that correspond to Ti and T2 respectively. 
For every such edge ct, the sets R{Ti) and R{T2) are disjoint regions, such that for i = 1,2, 
all vertices and edges of Gi are drawn entirely in the closure of Ri. 

The above lemma shows that T* is already close to a surface split decomposition of G' , 
and therefore of G. The fact that it is not actually a branch decomposition yet (since T* is 
not necessarily ternary and since non-leaves may correspond to edges of G) will be addressed 



later. First we give a bound on what will be the width of the resulting branch decomposition, 
in Lemma [HI For this we need the following simple fact. 

Proposition 13 Let g be the genus of the map G, vr. For and Ts as constructed above, it 
holds that \E{T+)\ - \E{Ts)\ = 2g. 

Proof: Denote x = \E{T^) \ — \E{Ts)\, and denote by n, m and / the number of vertices, edges 
and faces respectively of TZq'- By Euler's formula, 2 — 2g = n — m + f. (Recall that G, and 
therefore TZc can be 2-cell embedded in Sg.) The number of edges and arcs of is n — 1 + x, 
since x edges were added to a spanning tree of TZq'. The number of edges of T* is / — 1, since it 
is a spanning tree of the dual graph of TZg'- For every edge e of TZg', either e appears in Tg or 
its dual is part of T*. Therefore, m = (n — l+x) + {f — l). This yields 2 — 2g = n — m+f = 2 — x, 
and thus x = 2g. □ 

The following definitions are used in the proof of Lemma [TH An Euler tour is a walk that 
contains every edge of a graph exactly once. It is well-known that every connected graph in 
which all vertex degrees are even admits an Euler tour. A directed edge from vertex u to u is 
called an arc, and denoted by {u,v). An orientation of an undirected (multi-)graph is obtained 
by replacing every edge uv either by a arc {u, v), or an arc {v, u). In a partial orientation, only 
a subset of the edges is replaced by arcs. This results in a mixed graph, which contains both 
edges and arcs. The in-degree of a vertex v in a mixed graph is the number of arcs {u, v) (so 
incident edges are not counted). A diwalk of length in a mixed (multi-) graph G is a sequence 
vq, ei,vi, e2, ■ ■ ■ ,Vk where Vi € V{G), and e.j G E{G) is an arc of G from Vi-i to Vi. It is a 
dicycle if ^"0, . . . , Ufc-i are distinct, but vq = v^- 

Lemma 14 For an edge ex G E{T*), let Ti be the vertex set of one of the components of 
T* — ct, and let G'^ be the subgraph of G' that corresponds to Ti. For every such edge ex, there 
are at most (2(7 + l)(4(i + 3)/2 vertices of G'l that lie on the boundary of R{Ti). 

Proof: Consider the set J- of faces of TZc that corresponds to Ti . The boundary of the region 
R{Ti) is formed by vertices and edges of TZg' that are incident both with faces from and 
faces not in F. Denote by Vb C ^(T^g/) and Eb C E{TZg') these sets. Thus {Vb-,Eb) is a 
subgraph of TZg'-, which we denote by B. We give an upper bound for \Eb\, which will then 
yield the stated upper bound on the number of vertices in V{G') n Vb- 

Clearly all vertices in Vb have even degree in B. Therefore we can construct a collection 
Wi, . . . , Wc of Euler tours, one for every component of B. We show that the length of such an 
Euler tour Wi is bounded by x(4d + 3), where x is the number edges in that are not part 
of Ts- Edges of B that are not part of Tg will be called special edges- There are two kinds of 
special edges: firstly, there is a single edge e^ € E{1Zg') that is the primal edge of ex G E{T*). 
This is the only edge of B that is not part of . Hence all other special edges are edges in 
E(Tg)\E(Ts)- Therefore there are at most 2g-\-l special edges in total in B (Proposition llSp . 

Since Ts is a BFS tree of TZg> rooted at a vertex r € V{G') of eccentricity d (with respect 
to G'), the height of Ts is at most 2d + \: this holds because every edge uv in a walk in G' can 
be replaced by a walk of length two in TZg' from u to v- In addition, one extra edge may be 
necessary when considering the distance to a face vertex of TZg'- (That is, a vertex that is not 
part of G'.) 

Now we consider the following partial orientation of -B: we orient the T^-edges away from 
the root r. To be precise, we start with an orientation of Ts such that the root has in-degree 
and all other vertices have in-degree 1. All edges of B that are in Ts are oriented this way, and 
all other edges (special edges) remain undirected. We conclude that this partial orientation of 
B satisfies the following properties. 



• the maximum length of a diwalk in i? is 2d -|- 1. 



• B contains no dicycles. 



• B contains no vertices with in-degree greater than 1. 

It follows that every Euler tour Wi contains at least one special edge, otherwise there would 
be a dicycle or a vertex with in-degree at least 2 in B. Let Wi = vq, . . . , vi_i, vq, and assume 
w.l.o.g. that ^0^1 is a special edge. A segment of Wi is a maximal subwalk W' = Vx, Vx+i, ■ ■ ■ ,Vy 
such that VxVx+i is a special edge and W' contains no other special edges. Hence the edges 
of Wi can be partitioned into x segments, where x is the number of special edges in Wi. By 
the above properties, the length of a segment is bounded by 4d + 3: in the extremal case, W' 
consists of a special edge, a dipath of length 2d + 1 oriented in the opposite direction of W' , 
and a dipath of length 2d + 1 oriented in the same direction as W' . (Reversing the direction 
again would yield a vertex of in-degree at least 2.) As we observed above, the total number of 
special edges in B is bounded by 2g + 1. It follows that \Eb\ < {2g + l)(4(i -|- 3). 

Since B \s a. subgraph of IZc and TZc has a bipartition determined by the face vertices 
and the original vertices of TZq' (recall that original vertices are vertices in V{TZgi) n y(G')), 
every edge in Eb is incident with exactly one vertex of G' . Every vertex of G' that lies on the 
boundary of RiTi) is incident with at least two edges from Eb- It follows that the number of 
those vertices is at most \Eb\/2 < {2g + l){4d + 3)/2. □ 

At this point, T* is already close to a low-width surface split decomposition of G, as shown 
by the Lemmas 1121 and I14[ However, T* is not yet a branch decomposition since it is not yet 
necessarily a ternary tree, and it does not hold that a vertex of T* is a leaf if and only if it 
corresponds to an original edge of G. We will now make some simple changes to guarantee 
these properties. We will use the following observation on the maximum degrees of T* . Recall 
that original edges of G' are edges that are also in G. 

Proposition 15 A vertex Vd € V{T*) has degree at most 2 if it corresponds to an original edge 
of G' , and degree at most 3 otherwise. 

Proof: First we prove that every vertex of T* that corresponds to an original edge e = uv oi 
G' has degree at most 2. Because new edges have been added on either side of e, e is incident 
with two 2-faces in G' . Hence in the radial graph TZq' of G' , the face / that corresponds to e 
is incident with the vertices u, v and two vertices x and y, with d{x) = d{y) = 2. Therefore, 
the spanning tree Ts of TZc contains at least one edge that is both incident with x and with 
the face /. The same holds for y. It follows that the spanning tree T* of the dual graph of TZq' 
contains at most two edges incident with Vd- 

Similarly, it follows that every vertex of T* that corresponds to a new edge e in G' has 
degree at most 3: the new edge e is incident to at least one 2-face in G' . So at least one of the 
four edges of the face / of TZc that corresponds to e is included in Ts- The claim follows. □ 

Second stage of the construction Using Proposition [T5l we can modify T* as follows to 
obtain a new tree T and mapping fi from L{T) to E{G): Firstly, for every vertex Vd € V{T*) 
that corresponds to a new edge e € E{G'), with d{vd) < 2: delete Vd if d{vd) = 1, and suppress 
Vd if d{vd) = 2. Hence the remaining vertices of T* that correspond to new edges have degree 
3. Subsequently, for every vertex Vd € V{T*) that corresponds to an original edge e € E{G'): 
if d{yd) = 1, then set ^{vd) = e. If d{yd) = 2, then introduce a new vertex u and an edge uvd, 
and mark ^{u) = e. By Proposition [15] this concludes all possible cases. 

Proof of Theorem [TTJ We prove that the tuple (T, ^u) constructed above is a surface split 
decomposition of G, of width {2g + l)(4(i -|- 3)/2. From the operations above it immediately 
follows that T is a ternary tree and is a bijection from the leaves of T to the original edges 
of G, so (T, /i) is a branch decomposition of G. Consider an edge ex € E{T), and the two tree 



components Ti and T2 of T — e. For i = 1, 2, let Gi be the subgraph of G (not of G') induced 
by the edge set fi{L{T) n V{Ti)). 

First consider the case that one of the trees, say Ti, consists of a single vertex. Then Gi 
consists of a single edge ec- In this case it is trivial to define the disjoint regions Ri and R2 
such that ec lies in the closure of Ri and all other edges lie in the closure of i?2. Clearly 
|mid(er)| < 2 in this case. 

Now suppose both Ti and T2 are non-trivial. Then by the construction of T, ct is also an 
edge of T*, or bt is the result of suppressing a sequence of adjacent degree 2 vertices in T*, 
that all correspond to new edges of G' . In either case, we easily select a corresponding edge 
e* € E{T*) that separates the vertices of T* into components T* and T^, such that for i = 1, 2, 
fj,{L{T) n V{Ti)) is a subset of the edges of G' that corresponds to T* . Hence Lemma [T2] shows 
that the desired regions i?i and R2 exist (recall that G is a subgraph of G'). Furthermore, 
Lemma [T^ shows that at most {2g + l)(4(i + 3)/2 vertices of G lie on the boundary of these 
regions (recah that V{G') = V{G)). Hence |mid(er)| < {2g + l)(4d + 3)/2. 

Using the given combinatorial embedding vr of G, it can be shown using standard arguments 
and data structures that all graphs used during the construction of (T, /x) can be constructed 
in linear time (these are G', IZc Ts, T*, Tg and T itself), together with their combinatorial 
embeddings and with data structures that keep track of the bijections between the various 
vertex/face/edge sets. □ 



6 Summary of the algorithm for bounded genus graphs 

In this section we show how to combine the ingredients of the previous sections to obtain our 
main results; we give an algorithm for counting the number of subgraphs isomorphic to P in a 
graph G of bounded genus, and an algorithm for listing all of them. Without loss of generality, 
we may assume that G is connected. We first use a relatively standard argument related 
to decomposing embedded graphs into layers that admit low width branch decompositions 
(Lemma [T6]) . Eppstein used this argument for Subgraph Isomorphism in planar graphs [T^ , 
and in bounded genus graphs [15]. It was also used by Dorn [8]. 

Throughout this section we use the following notations. We choose an arbitrary vertex 
r S V{G), and construct a BFS tree Tbfs rooted at r. Let d be the eccentricity of r. For 
i = 0, . . . , d, define Li C V{G) to be the set of vertices at distance i of r, which are the vertices 
at height i in Tbfs- We call this set a layer. Let Gj = G[Li U Lj+i U . . . U Lj]. 

Lemma 16 Let G,tt be a map of genus at most g, with a vertex r of eccentricity d, and 
subgraphs Gj as defined above. For i, j with < i < j < d, a surface split decomposition of G\ 
of width at most {2g + l)(4j — 4i + 7)/2 can be constructed in time 0{n'), where n' = \ V{Gl)\. 

Proof: For every i, let the graph Hj be constructed from G| by adding one vertex u, and edges 
from V to every vertex in Lj. This graph has genus at most g as well, since it can be obtained 
from the graph G[Lq U Li U . . . U Lj] by 

• first deleting all edges that are incident with at least one vertex in Lq U . . . Lj_i and that 
are not in Tbfs, &rid 

• subsequently contracting all remaining edges with both end vertices in Lq U . . . U Lj_i. 

Observe that contracting a set of edges that does not contain a cycle can be done while main- 
taining an embedding in S^. Clearly the vertex v has eccentricity j — i + 1 in so by 
Theorem [m a surface split decomposition of of width at most {2g + l)(4(j — i + 1) -|- 3)/2 
can be constructed in linear time. This easily yields a surface split decomposition of Gl of 
width at most {2g + l)(4j - 4i + 7) /2. □ 



We now present a novel algorithm for counting the number of P-isomorphs in G, even for 
the case where P is disconnected. Let c be the number of components of the pattern graph P, 
which are denoted by P^, . . . , P'^. For S C {1, . . . , c}, denote by the subgraph of P induced 
by the components with labels in S. Formally, P^ = P[Ujg5l^(P*)]. For every Gj, we consider 
the following coloring, which uses the color set {!,... ,j — i + 1}: the vertices of layer are 
all colored with color x — i + 1. For < i < j < d with j — i < k and S Q {1, . . . , c}, we define 
DPC](S') to be the number of colorful subgraphs of Gj that are isomorphic to P^ . In other 
words, this is the number of P'^-isomorphs that use all layers of Gj . 

Proposition 17 Let g be a constant. Let G, vr be a simple map of genus at most g, let P be a 
graph on k vertices with c components, and let DPC- (S) and Gj be as defined above. For given 
i, j with < i < j < d and j — i < k, and given sei S C {1, . . . , c}, DPCf- (S) can be computed 
in time 2*^('') • 0{n'), where n' = |F(G^')|. 

Proof: By Lemma [T6l we can find a surface split decomposition of G' of width w € 0{g{j — i)) C 
0{k) in time 0{n') (recall that g is viewed as constant). The layer-based coloring for Gj uses 
q < j — i + 1 < k colors. Hence by Theorem [TOl we can compute the number of colorful 
subgraphs of that are isomorphic to P^ in time 2^^^^ ■ 0{n'). □ 

For j G {0, we define DPT''(S') to be the number of subgraphs of Gg that are 

isomorphic to P^ . Note that these are not required to be colorful. In particular, if S = 0, then 
we define DPT-' (5) = 1. To avoid the discussion of trivial cases, we simply define DPC^(5) = 
if i < 0, and DPT^ (5) = if j < 0. 

Lemma 18 For every j € {0, . . . , d} and 5" C {1, . . . , c}: 

DPT>{S) = DPTJ~\S) + Y,DPT>~'^-\Si) ■ i)PC^L,+i(52), 

where the summation is over all partitions of S into Si and S2 with S2 7^ 0, and all x € 
{l,...,k}. 

Proof: We prove that every subgraph G' of G that is isomorphic to P'^ is counted exactly once. 
Consider such a subgraph G' . If G' contains no vertices of layer Lj, then it counts towards 
DPT-'~^(S') as well, and therefore it is counted by the first term. In the summation term, since 
DPC^_^_l^^(S') counts only colorful graphs and x > 1, only subgraphs that contain vertices of 
layer Lj are considered. Hence in this case, G' is counted exactly once. 

Now suppose that G' contains a vertex of layer Lj . Then it is not a subgraph of Gq~ , so it 
is not counted by the first term. Let x be the maximum number such that G' contains vertices 
of all layers Lj,Lj-i, . . . , Lj-^+i- Since P^ has at most k vertices, we have x S {1, . . . , k}. Let 
S2 Q S he the indices of components of G' that are subgraphs of G-j_^_^_^, and let iSi C be 

the indices of components of G' that are subgraphs of G^~^~^. By choice of x, Lj^x contains 
no vertices of G', so every component of G' falls in one of the two categories, and therefore 
Si, S2 is a partition of S. Now we have deduced the unique combination of x, Si and ^2 such 
that G' consists of a colorful subgraph of G^_^_,_]^ isomorphic to P^'^, and a subgraph of G;j~^~^ 
that is isomorphic to P^^. Hence in this case, G' is also counted exactly once, namely by the 
corresponding term in the summation. □ 

Theorem 19 Let g be a constant. Let G, vr be a simple map on n vertices of genus at most g, 
and let P be a graph on k vertices with c components. In time S'^ • 2'-"^^'> ■ 0{n) C 2'^^'"'^ • 0{n), 
the number of subgraphs of G that are isomorphic to P can be computed. 



Proof: In the first stage of the algorithm, we compute DPC;^ (S) for every i, j with < i < j < d 
and j — i < k, and every S C {1, . . . , c}. There are 2^ choices of S, so computing DPC^ for 
one particular choice of i and j and all possible choices of S takes time 2'^ ■ 2'-'^^^ ■ 0{n'), where 
n' = \V{Gl)\ (Proposition I17p . Since every vertex of G appears in Gj for at most 0{k'^) choices 
of i and j, this yields a total complexity of 2^ • 2<^('=) • k^ ■ 0{n) C 2" ■ 2^^^^ ■ 0{n) for the first 
stage. 

In the second stage of the algorithm, we use the recursion from Lemma [18] for computing 
DPT-^(5) for every j e {0, . . . , d} and 5 C {1, . . . , c}. For a given j, it takes time 3^^ • 0{k) to 
compute this for every S: there are 3^^ partitions of {1, . . . , c} into Si, S2 and S (the complement 
of 5), and k choices of x, and therefore the total number of terms in all expressions (for given 
j) is bounded by 3*^ • k. Computing this for every j (in increasing order) therefore takes time 
3^ • 0{k) ■ 0{d) C 3^= • 0{k) ■ 0{n) C 2^^^^^ ■ 0{n). 

Finally, the algorithm returns the value of DPT'^({1, . . . , c}), which is the total number of 
subgraphs of Gq = G that are isomorphic to P. □ 

Now we show how to extend the above counting algorithm to an algorithm for listing all 
isomorphs. For a connected graph P on k vertices, Eppstein [T3] gave an algorithm for listing 
all P-isomorphs in a planar graph G on n vertices, in time 

20(fciogfc) 

^O(i) . 772, where m 

is the number of P- isomorphs. He asked whether this could be extended to the case where P is 
disconnected. We answer this question, by extending the result to the case where G may have 
bounded genus and P may be disconnected. Furthermore, we improve the complexity. The 
proof of our result follows the general method sketched by Eppstein [14^ Lemma 3] , which uses 
three stages. We will call these the counting stage, backtracking stage and generation stage. We 
apply this to our dynamic programming algorithm from Theorem [191 In addition, we fill in 
some essential details of the analysis, related to data structures, which is necessary to obtain 
the claimed complexity bound. 

Theorem 20 Let g be a constant. Let G, vr be a simple map on n vertices of genus at most 
g, and let P be a graph on k vertices with c components. In time 3'^ ■ 2^^^^ ■ 0{n) + mk^^^\ 
all subgraphs of G that are isomorphic to P can be generated, where m is the number of such 
subgraphs. 

Proof: The entire listing algorithm has three stages: a counting stage, a backtracking stage, and 
a generation stage. For the counting stage, we first run the counting algorithm described in the 
proof of Theorem [T9l Recall that this algorithm first computes the values DPC^ {S) for every 
S C {1, . . . , c} and < i < j < d with j — i < k using the algorithm from Section [3l and then 
computes the values DPT-' (5) for every S C {1, . . . , c} and j G {0, . . . ,d} using the recursion 
from Lemma [m 

For the backtracking stage, we compute for every DPT-'(5') and DPC^(S') whether they 
contribute to the total number of isomorphs, or contribute for short. This notion is defined 
recursively as follows, and denoted by 0/1 variables CT^{S) and CC^(5), respectively. 

• We initialize CT''(5) = for every j and S, and CC^(S) = for every i, j, S. 

• Recall that the total number of subgraphs of G that are isomorphic to P is DPT'^({1, . . . , c}) 
Therefore, we set CT'^({1, . . . , c}) = 1. 

• Every entry DPT-' (5"), in particular DPT'^dl, . . . ,c}), is computed by a summation of 
DPT^'^i(5) and a number of terms DPT-'"^~i(5i) • DPC^_^_^i(52) (Lemma [H]). For 
every entry DPT-' (5) that contributes to the total (i.e. that has CT^ (S) = 1), and every 
non-zero term DPT^-^-i(5i) • DPq_^^i(52) in its sum, we set CTJ~^'~^(5i) = 1 and 
CC^:_^+i(52) = 1. If DPTJ-i(5) / 0, then we set CT^-\S) = 1 as well. 



In the next stage we will generate lists of graphs for every j, S with CT^ {S) = 1, and every i, 
j, S with CC^(S') = 1. We have now set these values such that every graph that is generated 
is also part of at least one P-isomorph in G. It is easily seen that this backtracking stage has 
the same complexity as the counting algorithm, which is 3^^ • 0{k) ■ 0{n). Here we only consider 
the part of the counting algorithm where the entries DPT-'(S') are computed, which is called 
the 'second stage' in the proof of Theorem [191 

Now we consider the last stage, the generation stage. Suppose first that for every i, j and 
S, where DPC^(5) = x and 00^(5") = 1, we have a list LCj{S) that contains the x different 
colorful subgraphs of G'l that are isomorphic to . (How to obtain this list will be discussed 
later.) For the generation stage, we essentially run the counting algorithm again. However, 
every time that we consider an entry DPT-'(S') = y with CT-'(S') = 1, instead of recomputing 
DPT-'(S'), we construct a list LT^{S) that contains the y corresponding subgraphs of Gq. We 
construct the list LT^{S) corresponding to DPT-'(S') as follows. First, we take a copy of the 
list LT^~^{S). Secondly, for every non-zero term DPTJ-^-^Si) • DPC^._^+i(52) in its sum 
(Lemma llSp . we combine every graph in the list LT^~^~^{Si) with every graph in the list 
LCj_^_^_i{S2)- (These lists have been computed before since we have set CT-'~^~^(5i) = 1 and 
CC■J_^_^_l{S2) = 1 in the backtracking stage.) The resulting list of graphs is appended to the 
list LT^{S). By doing this for every DPT-'(S') with CT-'(S') = 1, this procedure ends with the 
construction of a list LT'^dl, . . . , c}), which is the desired list of all subgraphs of G that are 
isomorphic to P. 

However, to obtain the stated complexity bound, with an additive term mk^'^^^ instead 
of nmk^^^\ we have to be careful here. Consider for example the case that G contains m! 
subgraphs that are isomorphic to P, and are contained in layer Li (i.e. are subgraphs of G\Li\). 
Let C = {1, . . . , c}. Then the list LT^{C) has length m' . If we copy this list every time, when 
constructing the lists LT^(C) for j = 2,...,d, then we use VL{dm!) copy operations. Both 
m! and d may be much larger than k, so this cannot be bounded by the stated complexity 
20(fc) . 0(n) + mfc'^(^). On the other hand, if for the computation of LT^{S) we start with 
the existing list LT^~^{S), and append entries, then the original list LT^~^{S) is not available 
anymore for later computations. 

Therefore, we initialize the list LT^(S) in constant time with a pointer to the list LT^~^{S), 
and later interpret this as the list LT^~^[S) itself. If the list LT^~^[S) itself consists of a single 
pointer to another list LT^ (S), and contains no further entries, then instead we let the pointer 
for LT^{S) point to the hst LT^ (S) as well. This ensures that finding an actual entry from a 
list takes only constant time; there is at most one pointer lookup for every actual entry that 
we find. 

Similarly, for every term DPT^'^^^^(Si)-DPCj„^+i(52) where 5i = (and thus DPT^^'-'^^iSi) = 
1), the resulting list of graphs is the same as LGj_^_^^{S2)- In this case, we append a pointer 

to LGj_^_^^{S2) to the list LT^{S) in a similar way. 

Based on this implementation, we can now prove the desired complexity bound. For every 
j and S, every operation involved in constructing the list LT^{S) that does not consist of 
combining non-empty graphs from different list takes constant time. Therefore, the complexity 
of all of these operations can be bounded by 3^^ • 2'^^'') • 0{n) again (with the same argument 
as in the proof of Theorem I19p . Because of the backtracking stage, we only construct graphs 
that will be part of P-isomorphs of G. Therefore, every operation that combines (non-empty) 
graphs from a list LT^~^~^{Si) with graphs from a list iCj_^^^(52) can be attributed to the 
construction of a P-isomorph. Every such operation increases the number of components in 
the resulting graph (we take the disjoint union of non-empty graphs in with graphs 

in Hence for every P-isomorph there are at most c < k such operations, which all 

take time polynomial in k. Hence the complexity of all of these operations can be bounded by 
mk°^^\ This proves that the complexity of this part of the algorithm is bounded by 3*^ • 2*^^^^^ • 



0(?i) + mytO(i). 



It remains to consider the problem of computing a list of graphs corresponding to entries 
DPC^(S') with < j — i < k and 00^(5") = 1. That is, we have to construct a list of all colorful 
subgraphs of Gj that are isomorphic to P'^ . This can be done analogously, based on the counting 
algorithm from Section [3l The algorithm again consists of a counting stage, backtracking stage, 
and a generation stage. Let (T, fi) be the constructed surface split decomposition of Gl . For 
every e € E(T) and every entry (T, 7, A, rj) € Te, we first compute whether the entry contributes 
to the total number of colorful P'^-isomorphs. If so, we construct a list of r] subgraphs of Ge 
that correspond to this entry. Similar to the above argument, we use pointers. This way, every 
operation that results in a list of graphs that is the same as a previous list takes constant 
time, and can be attributed to the counting part of the complexity function. (This happens for 
instance when combining the empty graph from a list 7/ with graphs from Tg, where / and g 
are the children of e.) Every operation where graphs in lists are modified is attributed to one 
of the final colorful isomorphs. 

With a similar argument as before we can guarantee that only polynomially many operations 
are attributed to every such isomorph: When we combine entries such that modifying a list is 
necessary, in every case the number of vertices plus edges increases, or the number of 7-images 
decreases. This yields the stated complexity. □ 



7 Extensions and improvements 

Non-orientable surfaces Although we restricted ourselves to orientable surfaces for ease of 
presentation, we remark that our results hold for non-orientable surfaces as well. We shortly 
mention the facts about non-orientable surfaces that have to be used for adapting our proofs to 
this case. For c > 1, denote by Nc the surface obtained by adding c crosscaps to So (see [18] for 
details). Nc is called the non-orientable surface of genus c. For every non-orientable surface S, 
there exists exactly one value c > 1 such that S is homeomorphic to Nc \18\ p. 82]. For a graph 
G with n vertices and m edges that is 2-cell embedded in Nc with / faces, Euler's formula states 
that n — m + f = 2 — c. For graphs G that are embedded in Nc, the dual graph G* can be 
defined similarly as before, with one vertex for every face of G, and the same number of edges 
as G. (See [181 p. 103] for details.) Similarly, the radial graph TZg of G can be defined. With 
these ingredients, the constructions, bounds and proofs of Section [5] can easily be generalized 
to the non-orientable case. 

With regard to the bound in Section U Bender and Canfield [3] considered non-orientable 
maps as well, and proved that there are at most 0{m^^^~'^^^^12^) simple rooted maps of non- 
orientable genus c on m edges. We remark that two non-orientable embeddings of a graph G 
in Nc are called equivalent if facial walks are mapped to facial walks; we cannot define this 
based on clockwise edge orders now. Equivalence of non-orientable maps is defined similarly. 
Although Miller [17] treated only the orientable case explicitly, Grohe [16] analyzed a different 
isomorphism testing algorithm for both the orientable and non-orientable case, and proved that 
it decides isomorphism correctly in time n'^^^^ and n^^^\ respectively. This shows how the 
results from Section U] can be generalized. 

In Section [H the only observation related to embeddings was that contracting a set of edges 
that does not contain cycles can be done while maintaining an embedding in Sg. Clearly, this 
holds for Nc as well. In the remaining sections, no information about the embedding was used. 
We conclude that our results hold for the non-orientable case as well (with different constants). 

Induced subgraphs We observe that with a trivial modification, our algorithm can find and 
count induced subgraphs as well. Consider the dynamic programming algorithm from Section [3j 
For every edge e G E{T) of the branch decomposition, we now define the table Te to be induced 



k-correct if for every entry {H,^, A,ifj G there are exactly r] induced subgraphs G' of Ge 
with < k that are (i?, 7, 74)-subgraphs, and Te contains no pairs of equivalent entries. 

Note that it should be an induced subgraph with respect to Ge, not necessarily with respect to 
the entire graph G. For leaf edges e G E{T), this property is easily guaranteed: in this case, 
Ge is a graph that consists of a single edge xy. When initializing the table Te, we omit the 
subgraph that consists of two vertices and no edges, and include the other four subgraphs of 
Ge- This yields an induced /c-correct table. Observe that the dynamic programming update 
rule from Section [3] preserves induced fc-correctness, without further modifications. (However, 
we remark that a table may contain graphs that do not correspond to induced subgraphs of G; 
G may contain an edge between two vertices of Ge that is not part of Ge-) Hence Theorem [TOl 
holds for Induced Subgraph Isomorphism as well. The algorithm from Theorem [19] does not 
have to be modified either: it combines induced subgraphs of G'j_^_^_^ with induced subgraphs 

of G;j"^^^. Since these graphs are separated by the layer Lj^x, the result is again an induced 
subgraph. 

Planar graphs For bounding the table size in Section HI we applied a bound by Bender and 
Canfield [3], which is expressed in the number of edges of the graph. In the case of planar 
graphs, we can apply a result by Tutte, based on the number of vertices of the graph, to 
obtain better constants. Tutte [2^ showed that there are 0(9.5") equivalence classes of rooted 
triangulated maps on n vertices. We can modify the encoding in the proof of Lemma [7] as 
follows: given an entry {H,^, A,r}), we can embed H such that all 7-images are incident with 
the outer face, introduce the vertex u again, drawn in the outer face, and add edges from u to 
every 7-image in H. Choose one of these edges to be the root, similar to before. Next, we add 
edges until a triangulation H' , vr' is obtained. The edge labeling A is constructed by marking 
the added edges with 0, and the original edges with 1. This yields an encoding H' , vr', (u, v), A, p 
of the pair H,"f, where p and {u,v) are defined the same as before. Then H',7r' is a rooted 
triangulated map on at most fc + l vertices. Hence it has at most 3A; — 3 edges. Applying Tutte's 
result [24] and assuming that only one color is used (which may be done for the connected case, 
see [H]), it follows that there are at most 2"" ■ 0(9.5^) • 23^-3 = 2"^ • 0(76'') entries in the table. 
Although this constant is still very large, it is significantly better than the constants from the 
previous algorithm for planar graphs [8]. This is a strong indication that when implemented, 
our algorithm is likely to be much faster as well. We believe that for small values of k and 
connected graphs P, our algorithm may in fact be fast enough in practice, contrary to what our 
rough asymptotic complexity bounds suggest. Developing and testing a good implementation 
would be a good topic for further research. The key to a fast implementation would be a faster 
dynamic programming update step, based on ordered tables (see [13 [7] for examples). 
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